Non-local effects at the onset of the chiral transition 
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Inspired by analytic results obtained for a systematic expansion of the memory kernel in dissipative quan- 
tum mechanics, we propose a phenomenological procedure to incorporate non-markovian corrections to the 
Langevin dynamics of an order parameter in field theory systematically. In this note, we restrict our analysis to 
the onset of the evolution. As an example, we consider the process of phase conversion in the chiral transition. 
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I. INTRODUCTION 

For high enough values of temperature, strongly interacting 
matter should be in a quark phase due to asymptotic freedom. 
For the last few years, relativistic high-energy heavy ion col- 
lisions have been providing valuable information on the new 
state of matter that seems to have been created according to re- 
cent data from experiments at BNL-RHIC H], even though its 
true nature is still uncertain. Moreover, comparison between 
temperature estimates from data and results from Lattice QCD 
apparently shows that the critical temperature {Tc '--^ 150 MeV) 
for the chiral transition has been exceeded |2|. 

As the deconfined plasma created in a heavy-ion collision 
expands, it cools down, undergoing a phase transition (or a 
rapid crossover) back to hadronic matter, as well as break- 
ing the approximate chiral symmetry. Results from CERN- 
SPS and BNL-RHIC feature an explosive behavior in the 
hadronization process and seem to favor a very fast spinodal 
decomposition as the mechanism of phase conversion. Ef- 
fective chiral field theory models for QCD also point to an 
explosive spinodal decomposition scenario |3, 4]. 

Under these conditions, since relevant time scales are not so 
clearly separated (as, for instance, in the case of the primordial 
quark-hadron transition), the competition between several dy- 
namical phenomena will significantly affect the phase conver- 
sion process. In this context, fluctuations and non-local effects 
will play a major role in the dynamics of the order parameter. 
Therefore, a linear response non-equilibrium field theory ap- 
proach will provide a highly non-local Langevin equation for 
the evolution of the chiral condensate |5|, and the memory 
kernel details will presumably be important in the determina- 
tion of the physical time scales for the whole conversion pro- 
cess |a|. Most likely one will find interesting transient effects 
and, eventually, a hierarchy of time scales. 

Since the structure of memory integrals and colored noise 
that appear in field theory, even for relatively simple effec- 
tive models, is often rather complicated, one has to resort to 
numerical simulations to obtain exact results |7]. Analytic re- 



* leticia @ if . ufrj . br 
tfraga@if.ufrj.br 

* tkodama® if.ufrj .br 
^gkrein@ift.unesp.br 



suits, which could bring some qualitative understanding of the 
mechanisms involved, can only be achieved within systematic 
approximations 1 8J . 

In this paper, we propose a phenomenological procedure 
to incorporate non-markovian connections to the Langevin 
dynamics of an order parameter in field theory systemati- 
cally. We follow a path inspired by analytic results obtained 
for a systematic expansion of the memory kernel in dissipa- 
tive metastable quantum mechanics, a much simpler problem 
which we use as a toy model. To avoid some technical compli- 
cations, we restrict our analysis to the onset of the evolution, 
i.e., to very early times. Moreover, for simplicity, we inten- 
tionally neglect inhomogeneity effects and the noise contri- 
bution, although we know they will be important and should 
be incorporated IQl flflll . These issues, as well as a detailed 
comparison to exact (numerical) results, will be addressed in 
a future publication |8]. As an application, we consider the 
process of phase conversion in the chiral transition, compar- 
ing our results to those from a markovian Langevin evolution 
1 6] . In view of the caveats mentioned above, our results for the 
initial dynamics of the order parameter should be taken as a 
hint to the complex transient regime that can only be unveiled 
by a more rigorous approach. 

The paper is organized as follows. Section II presents the 
method in the framework of non-markovian dissipative quan- 
tum mechanics, giving special attention to the features that 
will be important in the phenomenological extension to field 
theory. The early-time dynamics of the chiral condensate with 
the first phenomenological non-local corrections is discussed 
in Section III. Section V contains our conclusions and outlook. 



II. NON-MARKOVIAN DISSIPATIVE METASTABLE 
QUANTUM MECHANICS 

The development of a systematic procedure to take non- 
local effects into account should have as starting point a care- 
ful determination of the domain of validity of the local approx- 
imation. However, as pointed out above, the structure of mem- 
ory kernels that come about in non-equilibrium field-theoretic 
approaches is rather intricate. This may turn the determination 
of the domain of validity into an excessively complex task. To 
isolate the analysis of the role of non-locality, we consider a 
much simpler case in quantum mechanics (QM), where the 
problem can be solved in a detailed and controlled fashion. 



Our QM toy model is constituted by a particle, with co- 
ordinate q and momentum p, interacting linearly with a heat 
bath which introduces dissipation in the dynamical evolution. 
The reservoir is modeled in the usual Caldeira-Leggett fashion 
lllll . The hamiltonian is thus given by: 



more and more non-local corrections. Taking advantadge of 
the localized profile of the kernel that appears in Eq. (0} and 
assuming that the system will eventually thermalize in a stable 
vacuum, we can derive the following series for the memory 
integral (for details, see [SJ): 



H = ^+V{q)+HR + Q{t)qY^CkRk; 
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where the sum over the heat bath modes is limited by a "cut- 
off" n, that corresponds to the maximum frequency for the 
interaction of the system with the medium. We also assume 
the system to be suddenly immersed into the reservoir at f = 
(which represents a rough simulation of the system "quench" 
in the chiral transition), so that the approximation of an ini- 
tially uncorrelated density matrix is reasonable: 



P(f<o) = (|o)(o|X,® 
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where T = 1/p is the reservoir temperature. 

The evolution towards equilibrium is studied within the 
Schwinger-Keldysh closed-time path framework |12|. Inte- 
grating over the bath variables, we obtain the effective action. 
After performing a Wigner transform to separate the classi- 
cal coordinate, Q{t), from the fluctuations (associated with a 
Gaussian noise ^(f)) the equation of motion assumes the fol- 
lowing semiclassical generalized Langevin form: 
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where r| is a phenomenological dissipation parameter encod- 
ing information about the heat bath and V (2) is the derivative 
with respect to Q of the modified particle potential. 



v'{Q)=V'{Q)-[2^a/K]Q. 



(5) 



In the markovian limit, this correction represents an ultra- 
violet divergence that is associated with the unphysical situ- 
ation in which the particle interacts with all frequencies of the 
heat bath. Usually, in the context of quantum dissipative sys- 
tems, this alteration in the potential is ignored, under the ar- 
gument that the interest resides on the modifications brought 
about by dissipation and not by medium renormalization ef- 
fects. In quantum field theory, however, the "particle" and 
the reservoir are deeply entangled, corresponding to different 
Fourier modes of the same field (refer to j^l for further discus- 
sion). Since we are ultimately interested in the chiral transi- 
tion case, we compare in this analysis the evolutions with and 
without the modification in the potential. 

The limit in which Q. is strictly infinite yields a local dis- 
sipation term in Eq. 0. Therefore, the non-markovian con- 
tributions arise from the assumption that Q. is still large, but 
finite. As one goes down in Q. one should, in principle, include 
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where the time-dependent coefficients J„{Q.t) contain infor- 
mation about the specific form of the original kernel. It should 
be noticed that the first term corresponds to the usual marko- 
vian approximation, in the limit i2 — > oo. 

The oscillatory behavior of the integrand in Eq. © is inher- 
ited by the J„'s. Nevertheless, Jo and Ji have definite signs: 
7o > and Ji < 0. However, for n > 3, one can show that 
the signs of those coefficients change according to the time at 
which it is being evaluated. Even though the amplitudes of 
the associated terms decrease, rendering these oscillations ir- 
relevant for long times, the transient dynamics will be greatly 
affected. 

We remark that the localized profile of the memory ker- 
nel is the essential feature for the development of the non- 
local series described above. The oscillatory characteristic of 
the 7,,'s is a specificity of the function sin(.Ji:)/ji: and its con- 
sequences do not encode modifications due to non-locality. 
Therefore, the alternating signs should not be interpreted as 
intrinsic memory effects. As a matter of fact, the essential 
feature being added here, as compared to the markovian case, 
is a non-zero correlation for times separated by an interval 
Af < 1 /i2. For ii large enough, higher-order contributions in 
Eq. (|6j become more and more negligible and the series can 
be truncated. This procedure systematically incorporates the 
effects from correlations between times within this interval. 

Now we analyze the non-local Langevin evolution dictated 
by Eq. (|6|l for the double-well potential, namely for V'{Q) = 
4A-2^ —MkiqQ, at early times. In this regime, the equation is 
linear For simplicity, we also approximate the coefficients by 
constants. Keeping terms up to the third-order time derivative, 
one has: 



MQ - [M&l + 2rin/jt] G + ri e + ri3 e= ^ (f ) (7) 
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The solution of Eq. Q is obtained through a straightfor- 
ward calculation, using the Laplace transform. Results for the 
average {Q) as a function of time are displayed in Fig. ^ for 
(2(0)} = (2(0)} = 0. One can see that both memory and po- 
tential corrections tend to accelerate the initial decay process 
towards the minimum at 2 = Qvac- It is clear that the modifi- 
cation of the potential introduces the major effect, remarkably 
reducing the time scale involved in the decay process. Al- 
though this correction may be ignored in the analysis of dissi- 
pative phenomena in QM, as discussed previously, it must be 
taken into account in the field theory case. In the latter, the 
medium renormalization of the effective potential should play 
a significant role in the phase transition process. 
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FIG. 1 : The average of Q in units of its vacuum value as a function of 
time for M = coq = "n = 2 in arbitrary units and (2(0)) = 0.04 Qmc- 
The dashed curve corresponds to the non-markovian case with fl = 
2 A ~ 3.852, while the dot-dashed one is the non-local solution in- 
cluding the potential renormalization for the same value of fi. 



The fact that not only the medium modification of the po- 
tential but also the memory correction tend to speed up the 
decay of the particle seems surprising at first sight. In fact, 
one expects the latter to slow down the rollover process. The 
behavior observed in Fig. ^is certainly a transient effect and 
should be washed out in the long run by dissipation and non- 
linear effects 1 8]. Nevertheless, it points to the possibility of 
a non-trivial transient regime at early times. We will return to 
this point in the next section. 



III. CHIRAL TRANSITION EARLY-TIME DYNAMICS 

Inspired by the results obtained above for the QM toy 
model, we can construct a phenomenological generalization 
to the case of field theory. In particular, we consider the case 
of the chiral transition. The effective potential, Ueff{^^T), 
for the chiral condensate i^{x. t) can be obtained, for instance, 
from an effective chiral model as described in Ref. [4]. 

The local form of the Langevin dynamics for the chiral con- 
densate, which plays the role of a non-conserved order pa- 
rameter in a Landau-Ginzburg approach, can be derived un- 
der convenient approximations within the framework of non- 
equilibrium field theory |5|l, and is given by 



n^+r(r)^ 



di 



-Ulff{ifJ)^^{x,t). 



(9) 



The second term in Eq. (|5Jl corresponds to markovian dis- 
sipation and the right-hand side represents a gaussian, white 
stochastic force, both attained, in principle, after drastic ap- 
proximations of complicated memory kernels 1 5]. Some of us 
have recently analyzed the role of dissipation in the hadroniza- 
tion of the quark-gluon plasma under these conditions |6]. 

The inclusion of non-local effects is done, phenomelogi- 
cally, by inserting higher-order time derivatives of (|)(x,f ) into 
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U',fA^J)^i,{x,t). (10) 



Trying to follow the scale hierarchy that arises naturally in 
the series expansion in QM, we choose the form of the coef- 
ficients accompanying these new terms. In addition to being 
proportional to the phenomenological temperature-dependent 
dissipation coefficient F, which can be estimated using kinetic 
theory, the F„'s are associated with increasing inverse pow- 
ers of a large scale Q.: T„{T) = (-l)«F/n"-^ Here, Q. is 
the parameter that controls the revelance of non-locality and 
imposes the intended hierarchy between the contributions of 
different orders of time derivatives. For the above truncation 
of the series that presumably appears in Eq. jlO> to be valid, 
i2 must be large (although finite) compared to the temperature 
r, the natural scale of the system. 

To obtain analytical solutions of Eq. (llOt . one can resort 
to a sixth-order polynomial fit of Ueff{(^,T) fff\ and investi- 
gate the short-time dynamics. Since we consider the explosive 
spinodal decomposition scenario, the potential is assumed to 
be time-independent. Its form is fixed at the spinodal tem- 
perature (Tsp ^ 108 MeV), and the system is initially slightly 
displaced from the inflexion point ((|)o ^ 0.162 Tsp). 

For short times, f/^^^ can be approximated by its linear form 
around (|)o. We treat, for simplicity, the dynamics of the mean 
value of the field over the stochastic noise for a homogeneous 
initial condition. Therefore, Eq. (I10> becomes equivalent to 
an infinite set of linear ordinary differential equations: each 
Fourier mode of the field evolves according to an analogue of 
Eq. 0. When one examines the solution averaged over space, 
one verifies that the only relevant mode is k ~Q. Under these 
conditions, the field theory case is reduced to a QM problem 
for the Fourier component (^k=Q = ^p: 
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In Fig. 121 we present results for the initial dynamics of the 
order parameter (|)(x,f ) subject to a non-markovian Langevin 
evolution. The behavior at times larger than those shown in 
the figure reproduces qualitatively the QM case. Fig. ^ 

The existence of a competition between effects from dif- 
ferent memory coiTections is illustrated in Fig. |2] Although 
these results are valid only for short times, the study of these 
curves enlightens the role of each term in the dynamics of the 
phase conversion. 

The early influence of non-local dissipation is to retard the 
phase transition process. At extremely short times, the lower- 
order time derivatives are subdominant, so that (p w e with 
A- ~ {b/aY'^ o= ii^'^. Thus, the outcome of smaller values 
of Q. will be a greater delay in the initial evolution. On the 
other hand, for intermediate times, the modification brought 
about by the third-order derivative is overcome by the non- 
local correction of the potential. In this case, as Q. increases, 
the exponent X will decrease, X w /?'' 2, and the curves tend to 
the markovian one. 
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FIG. 2: The evolution at T = T^, of the average value of the order pa- 
rameter over noise and space in units of its vacuum value. In addition 
to the local evolution (full line), three curves with non-markovian 
corrections are shown: Q./T = 10, 5 and 2. 



Features that were neglected in this study should play a 
major role in the transition dynamics, possibly washing out 
those transient effects observed here. We will discuss the main 
shortcomings of the method in our conclusion. 



IV. CONCLUSIONS AND OUTLOOK 

We constructed a phenomenological generalization of the 
Langevin dynamics experienced by an order parameter in field 
theory, accounting for non-markovian correlations in a sys- 
tematic fashion. We used as inspiration a toy model in dis- 
sipative metastable QM, where the domain of validity of the 
local approximation was identified, allowing us to derive a 
derivative expansion of the memory kernel. In this series, the 
first term corresponds to the usual markovian assumption and 
higher-order time derivatives are accompanied by increasing 
powers of the scale 1 /ii, associated with the maximum time 
interval in which correlations exist. 

It should be noticed that this systematic expansion method 
developed can not be applied to any memory kernel. In fact. 



the kernel must have a localized profile, so that correlations 
between sufficiently distant times vanish. This is in general 
not true when one considers a semiclassical expansion around 
a non-trivial background as, for instance, in the case of dy- 
namical viscosity in nucleation 1 13|. 

Results presented for the onset of the non-markovian evo- 
lution of the chiral condensate neglecting stochastic noise ef- 
fects revealed an intricate transient structure caused by the 
competition between different contributions. The influence of 
non-locality at extremely short times is to delay the transtion 
process through a dissipation term proportional to a higher- 
order derivative of the field. At intermediate times, non- 
markovian modifications of the potential dominate the evo- 
lution, accelerating it. 

In this work, we identified the role of each non-local coiTec- 
tion in a very simple model for the early dynamics. The com- 
plete evolution will be studied numerically in a future publica- 
tion 1 8]. In a more realistic approach, however, several other 
phenomena should be taken into account and an even more 
complex balancing will dictate the evolution towards the new 
vacuum state. The inclusion of the complete non-linear po- 
tential, for instance, would accelerate the early evolution and 
slow down the process for larger times, smoothening or even 
erasing the initial modifications observed above. The effect 
of stochastic noise should be considered, as well as inhomo- 
geneities in the effective potential. Not only temporal, but also 
spatial non-locality should be taken into account. All these as- 
pects could, in principle, affect significantly the time scales of 
the phase conversion. 

To obtain reliable quantitative results in the field theory 
case, one must firstly derive this phenomenological non-local 
Langevin equation through non-equilibrium quantum field 
theory. In this procedure, the medium renormalization of the 
potential (which provides significant corrections in the QM 
case) would be incorporated naturally and the scale Q. could 
be identified and, possibly, estimated. 
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